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1. Introduction 

Models of statistical mechanics defined over Bethe lattices [I] constitute a framework 
that, due to two peculiar ingredients, namely exact solvability and finite connectivity, 
as opposed to unsolvability of finite dimensional models (for D > 2) or to exact 
solvability of models defined over the fully connected graph (mean field), sheds some 
light toward the understanding of more realistic models. It is in fact believed that 
several thermal properties of the model could persist for regular lattices as well. 
Nowadays, the emerging science of networks has attracted a renewed interested in 
such models [31 |31 [H H]. In fact, roughly speaking one can say that the Bethe 
lattice represents the simplest prototype of network where an exact solution is often 
accessible. At least for uncorrelatcd random networks, it has been indeed proved 
that, if we indicate with < f{q) > the average over the network of some function / 
of the vertex degree q, the critical behavior of a given model with positive couplings 
defined over a Bethe lattice of degree qs, coincides with that defined over the network 
provided that an effective substitution qs ^< q^ > / < q > tor taking into account 
of the distribution of link ends rather than that of links, is performed [SlITj. 

As is known, the exact solvability of a model defined over a Bethe lattice relies on 
the fact that the lattice is a tree, i.e., a graph without loops [1:. However, in relatively 
recent years progresses have been made in facing, analytically, models defined over 
lattices in which a finite number of closed paths per vertex is also allowed , which 
we refer as "generalized tree-like structures" , sec Fig. [T] and models defined over 
graphs in which the tree-like structure is partially broken due to the presence of an 
infinite number of closed paths per vertex as happens in generalized Bethe lattices (see 
for example [TU] and references therein) and in Husimi trees [TTl [TH [131 [HI [121 [H] . 
The graph depicted in Fig. [2] is an example of a Husimi tree. 

Notice that, by definition, in a loop both vertices and bonds repetitions are 
forbidden, whereas in a closed path only bonds repetitions are forbidden; in a closed 
path a same vertex can be crossed more than one time. So that, for example, in the 
graph of Fig. [l] and in a regular /^-dimensional lattice with D > 1 or in the graph 
of Fig. [21 the number of loops per vertex is finite, but the number of closed paths 
per vertex is finite only in the graph of Fig. [T] This kind of difference reflects on 
the different complexities in solving models defined on generalized tree-like structures 
and non tree- like structures; whereas the first class of problems can be exactly faced 
in the framework of a suitable Bethe-Peierls approach (see for example the reviews 
[151 [TOl ) or, alternatively, by using the rules given in the Refs. [SI [S] for determining 
the critical surfaces, for the second class of problems, in general, a closed analytical 
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Figure 2. The graph over which the competitive model of Eq. [T]is defined: a 
triangular Husimi tree. The continuous segments constitute a Bethe lattice F (an 
infinite tree). Dashed segments connect vertices of F belonging to the same shell. 
Note that, a part the root 0, the degree of F is = 3, whereas the degree of the 
Husimi tree is q'j, = 4. 

method is lacking. An exception to this is provided by models defined over generalized 
Bethe lattices and Husimi trees. For these models in fact, the problem of statistical 
mechanics can be mapped on a that of a dynamical system involving few degrees of 
freedom [TUIIII]. In particular, the model defined over the Husimi tree of Fig. 
which is the subject of this paper, can be solved analytically in a closed form. 

In these kind of relatively solvable models, one is usually interested in studying 
competition effects, so that two or more independent couplings Ji,J2,... along 
the links of the lattice are given. As already investigated by numerical and 
semi-analytical methods, not only on generalized Bethe lattices and Husimi trees 
[HI HOI ini HI US Eni m, but even on finite dimensional lattices [13 [5^, due to 
the competitiveness in amplitude and sign of the couplings and to the presence of 
loops giving rise to frustration [27], a quite rich scenario of phases can take place. 
In fact, besides ferro and antiferromagnetism, competing interactions lie at the heart 
of a variety of original phenomena in magnetic systems, including the existence of 
modulated [TU] and spin glass phases US] . 

After the study of the Ising and the Potts competitive models O [H] [HI [TT] 
defined over generalized Bethe lattices and Husimi trees with a uniform or periodic 
distribution of the couplings Ji , J2 , . . . , the next natural step toward the complexity 
of these competitive models is the introduction of disorder. We present here a twofold 
study of the simpler Ising model with two competitive interactions Ji and J2 defined 
over the triangular Husimi tree of Fig. [2] We first review this model giving the 
proper interpretation of the phases in terms of the magnetization. Then we introduce 
the disorder to study the triangular ±J model. To this aim, we let Ji and J2 to 
be random couplings having a discrete probability distribution whose spread, i.e., 
the measures of the disorder, is characterized by two parameters, pi and p2- The 
analysis of the phase diagram of the random model is done by using an effective 
method recently developed in [55] and already used in [30] ■ This method consists on 
a mapping between the phase diagram of the random model and the corresponding 
non random one. In |29| it has been shown that, in the thermodynamic limit, the 
upper phase boundary of the random model one finds by using this mapping, as well 
as the corresponding critical behavior, become exact whenever the graph over which 
the model is defined turns out to be infinite dimensional in a broad sense, as happens 
for instance in generalized tree-like structures or generalized Bethe lattices and Husimi 
trees. Unlike what emphasized in a previous version of this manuscript however, we 
stress here that Bethe lattices (here meant as infinite trees) and Husimi trees represent 
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a limit case where the mapping in general is not well defined. In fact, the mapping 
can establish a correspondence between the random model and a corresponding non 
random one, both built up over the same graph, only if their density free energies 
exist (more precisely the non trivial part of the density free energy, see Eq. (|12l) '). We 
refer the reader to the recent work [31| for a rigorous formulation of the mapping. It 
happens unfortunately that, due to the problem of the boundary conditions (see also 
the discussion in Appendix) , which in a Bethc lattice and in a Husimi tree constitute 
a finite fraction of the total size of the system, in a Bcthe lattice (here meant as an 
infinite tree) and in a Husimi tree, the density free energy, as a proper thermodynamic 
limit, does not exist (however the circumstances for a Bethe lattice is in some sense 
fortunate and the mapping turns out exact). Therefore we advise the reader that 
the results we report in this manuscript cannot be taken as exact; they have to be 
considered as a mere application of the mapping being here formal and not rigorous. 
Nevertheless, the equations of the mapping are still defined and, as an approximation, 
they can be applied in a Husimi tree case as well. In fact, in studying another random 
version of the model, the diluted frustrated model, in which the disorder is introduced 
by randomly deleting the negative couplings J2 with some probability we have 
found some differences between the results one can obtain by applying the mapping 
and the results one can find by using a Bcthe-Peierls-like approach [32j . It turns out 
that, in the limit < — J2 ^ Ji, this latter method succeeds in giving the exact 
percolation threshold Hp = 1/^/2, while the mapping gives fip = VS— 1; and a larger 
difference is observed in the opposite limit — J2 ^ Ji > 0. Nevertheless we think that 
the application of the mapping also in the Husimi tree cases provides easily important 
information of the random model. We stress in fact that in the mapping no ad hoc 
approximation, like the local tree-like approximation, is introduced. 

For the present model, Eqs. UJH]), we have found the following scenario. 
Concerning the magnetization, in the non random model in the frustrated regions 
there is no phase transition for IJ2I > |Ji| even at T = 0, and this is physically 
explained by looking directly at the ground state that becomes equivalent to a gas of 
non interacting dimers with coupling J2 (also called a spin liquid in the literature). 
In the ±J triangular model instead, even for pi = p2 = I, there exists always a stable 
spin glass transition and, for suitable values of pi and p2 a ferro or antiferromagnetic 
disordered transition may exist as well. Whereas the latter transition can be read just 
as due to an effective renormalization of the couplings Ji and J2 , the former transition 
is interpreted as due to the onset of the overlap ordering of the dimers, as signaled by 
the study of the exact ground states of the non random model which is equivalent to 
the random one for pi = P2 = 1 • Furthermore, we observe that in the limit | J2 1 2> \Ji\ 
the picture in terms of non interacting dimers is recovered also in the glassy phase. 

We shall limit ourself mainly to the study of the critical surfaces inferring little on 
the nature of the several phases. Furthermore, wc recall that for negative couplings, 
the correspondence between the solution of the Ising model on the Bethe lattice (the 
graph of Fig. [2] with J2 = 0) and that on a corresponding sparse system (like the 
Erdos Rcnyi random graph), where long global frustrated loops are also present, is 
lost. Therefore, following the same philosophy of, e.g., the Refs. [201 HH HO], the 
solution we will present here concerns, strictly speaking, only an Ising model model 
defined over a structure where global loops are absent. Nevertheless, wc believe that 
since we take into account the effect of all the other explicit loops, at least for certain 
aspects, the solution wc find sheds some light also on a possible sparse representation 
of the model even in the frustrated regions. In fact, the ground state we describe in 
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Sec. 3 has a quite strong similarity to that one finds in the anisotropic triangular Ising 
model p4] . 

The competitive non random model was partially already analyzed in [13|, 1151 132j . 
Notice that it is different from the model considered in the Ref. [TT] (third item); in 
fact, in the non random case, its relative simplicity allows an exact analytical study. 
Note also that in the Refs. [T3J [13 [32| no detail about the ground states was given. 

Very recently, a renewed attention has been attracted by the role of loops both in 
statistical mechanics and computer science [35l [36l |37l |38]. In fact, even if the most 
used Bethe-Peierls approximation usually works well in many important models near 
to be tree-like, rigorously speaking, it becomes exact only in perfect tree-structures 
and generalized tree-like structures, like the one depicted in Fig. [H whereas in more 
complex graphs like that of Fig. [2 can be wrong. Successful efforts have been made 
toward a general treatment of the problem and systematic corrections to the Bethe- 
Peierls approximation have been now rigorously established [32] ■ However, as a general 
rule, we recall that between the complexity, both analytical or numerical, in solving 
a random model and a corresponding non random one there is a gap. In fact, in 
general the approach used to solve, possibly exactly, the latter, does not lend itself 
to be simply adapted to the former. Our study, even if limited to the knowledge of 
the upper critical surface and to the critical behavior, constitutes an example of a 
general procedure which, at least in the cases where the conditions for the mapping 
are satisfied, covers this gap. 

The paper is organized as follows. In Sec. 2 we introduce the uniform and 
the random model defined over the triangular Husimi tree of Fig. [H In Sec. 3 we 
analyze the phase diagram of the uniform model giving the proper interpretation of 
the invariant ground states. In Sec. 4 wc recall the mapping through which the phase 
diagram of the random model will be analyzed in the following Sec. 5. Finally, some 
conclusions arc drawn in Sec. 4. The Appendix A is devoted to Sec. 3. 

2. Competitive Ising models on a structure with loops 

Let us consider a semi-infinite Bethe lattice F in which any vertex, but a singular one 
0, has coordination number = 'i {i.e. branching number fc = 2) as in Fig. [21 Chosen 
as a root vertex, it is convenient to divide F in shells 5*1, 5*2, . . . , iSn, . . . [l]. Given 
a couple of vertices [x, y) of F, we will write {x, y) as < x,y > if x and y are nearest 
neighbors, i.e., if they are connected through a bond of F, and (x, y) as > x,y < if x 
and y belong to the same shall Sn and are at distance d(x, y) = 2 on F. 

2.1. Non random case 

Let us consider the following Ising model with competing interactions 

H{a) = - Ji ^ <j{x)a{y) - J2 ''(^Mv) (1) 

<x,y> >x,y< 

where Ji , J2 S R are coupling constants and cr is a configuration of spin variables 
= ±1 over the vertices of F. Far from the root the graph over which this Ising 
model is defined, is characterized by the fact that at any vertex is present a triangular 
loop with two couplings Ji and one coupling J2 as indicated in Fig. [21 For J2 7^ the 
model is defined over the Husimi tree of Fig. [H having degree q'j^ ~ 4. 
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2.2. Random case (spin glass) 

In the random version of the competitive model of Eq. ([T]), the couphngs Ji and J2 
become independent random variables so that the corresponding Hamiltonian for the 
quenched system is 

^M = - E J^H^^iv)- E 4X^)'^(y)- (2) 

<x,y> >x,y< 

The free energy F is defined by 

-PF^jdV {{Ji]l}, Wl}) log [Z ({JW}, Wl})) , (3) 

where Z {Jx^l}j is the partition fmiction of the quenched system and 

dV { Jify}^ is a product measure given by 

dv[{Ji]l},{Ji%}) = n d^^l{J^]l) n dfi2{4'l), (4) 

<x,y> >x,y< 

dfii and dn2 being two given normalized measures featuring the disorder. In general 
this model, besides a disordered ferromagnetic or antiferromagnetic transition (P- 
F/AF), may manifest a spin glass transition (P-SG). A generic inverse critical 
temperature of the random model will be indicated by Pc- We will consider the 
following choice: 



2.2.1. Triangular ± J model 

dfi^(^4]lj - Ji) + (l~Pi)5(jW + Ji), (5) 

dfl2 [4^.1) = P2<5(4'^ - J2) + (1 - P2)<5(4'^ + J2), (6) 

where Ji and J2, as in the non random model, are arbitrary parameters and 
S [0,1] characterize the disorder, maximum at pi = P2 = 0.5 and minimum 
for pi and p2 equal to or 1. 



3. Analysis of the non random model: phase diagram - gas of dimers 

Concerning the phase transitions, the model defined in Sec. 2.1, even in the presence of 
an external magnetic field h, has been recently solved analytically in |13| . In Appendix 
A we report another approach for deriving the phase transitions. There are two critical 
lines given by 



26*1 
26*1 



' 1 - 361 

where the parameters 9i and 62 are defined as 



h > V3, (7) 



9i < I/V3, (8) 



^e^f'^ 02 = 6^"'^ (9) 



In Fig. Owe show the corresponding phase diagrams obtained as solutions of Eqs. ([7]) 
and dS]) along the axis T/Ji and J2/J1. Equations ([7]) and dH) describe, respectively. 
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a paramagnetic-ferromagnetic (P-F) transition, in which Ji > 0, and a paramagnetic- 
antiferromagnetic (P-AF) transition in which Ji < 0. By a direct inspection of the 
minimum of TJ, Eq. ([1]), it follows that, in the two non paramagnetic frustrated (Fr) 
regions Ji > 0, J2 < and Ji < 0, J2 < 0, sec Fig. [31 the corresponding ground 
states, despite the frustration, are simply F and AF, respectively. The reason for that 
lies on the smallness (in modulo) of the coupling J2 in the two mentioned regions. 
More precisely, by looking at the configurations minimizing the energy, it is easy to 
see that, if e.g. Ji > and J2 < 0, for Ji > | J2I there exist only two ground states 
obtained by taking all the spins parallel, the corresponding ground state energy per 
triangle being e = — 2Ji — J2 (see Fig. HJ; whereas for Ji < | J2I the ground state 
becomes infinitely degenerate, the corresponding ground state energy per triangle 
being e' = J2 (see Fig. [5]). In the latter case the whole set of ground states is 
obtained by taking in all the possible ways the spins connected by the couplings J2 
(i.e., the spins on the same shell at distance d = 2 on F) as antiparralcl. A similar 
conclusion holds also in the other frustrated region where Ji, J2 < 0. In this case, 
when I Ji I > I J2 1 , the ground state is obtained by taking alternated spins on the Bethe 
lattice F (note that, as a consequence, the spins are not alternated on the Husimi tree; 
spins on the same shall are parallel, see Fig. [6]); whereas for |Ji| < IJ2I the ground 
state has the same structure of the case Ji > and | Ji| < | J2I above analyzed (Fig. 
[5|). In other words, for T = in the Fr regions, crossing the critical P-F or P-AF 
lines amounts to move oneself between a connected set of triangles with ordered or 
antiordered spins (on F), respectively, and a set of independent dimers having at their 
edges two antiparallcl spins interacting with a negative coupling J2. In the sectors 
II and IV of Fig. [3l at low temperature, the P phase corresponds therefore to a gas 
of non interacting dimers having each one two possible degenerate states. Note that, 
as a consequence, in the regions IJ2I > |Ji|, for any selected ground state a, the 
spatial average magnetization is zero, but the Edwards Anderson order parameter 
Qea = Qa.a is 1, so that, in certain aspects, the zero temperature limit in these regions 
leads to a sort of glassy phase at zero temperature [43]. We note that, in this phase, 
the interesting "objects" are not the spins, but the dimers. In fact, by looking at the 
moments of P{q), the probability distribution of the overlaps between spins [28| . and 
using the symmetry among the ground states one has P{q) ~ S{q), whereas for the 
overlap between dimers one has P^{q^) = S{q^ — 1), where for overlap between dimers 
q^ we mean qj^^ = l/NJ2i < ^i^l >a< <^i<^i >i3, where a and f3, with a ^ P, label 
two ground states, the index i stands for a dimer's index and the superscripts x and 
y label the two spins of a dimer. 

In the literature this phase was called a spin liquid phase |32| . An important 
question to be addressed is whether this zero temperature phase signals the existence 
of a spin glass transition at some finite temperature or it is merely a feature of the 
zero temperature limit. In the next sections by applying the mapping we will see that 
the former hypothesis emerges naturally. 

4. Mapping 

4.1. Equations of the mapping 



Let 



ti = tanh(/3Ji), 
t2 = tanh(/3J2). 



(10) 
(11) 
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T/Ji 




Figure 3. Phase diagram of tlie competitive model, Eq. The transition 

lines P-F and P-AF come from the solution of Eqs. JT)! and lO, respectively. Fr-F 
and Fr-AF stand for the two frustrated ferro and antiferromagnetic regions with 
Jl > and J2 < 0, and Ji < and J2 < 0, respectively. At T=0, the P-F and 
P-AF lines cross the J2/J1 axis at J2/JI = —1 and J2/J1 = 1, respectively. 




Figure 4. One of the two ground states of the competitive model, Eq. in 
the Fr-F region (Ji > 0, J2 < and Ji > | J2I). 




Figure 5. One of the infinitely many degenerate ground states of the competitive 
model, Eq. ([TJ, in the regions Ji > 0, J2 < or Ji < 0, J2 < with | Ji| > | J2I. 
This phase can be read as a gas of non interacting dimers with antiparallel spins. 
A dimer here is defined as constituted by two spins belonging to the same shall 
and connected through a bond with coupling J2. 
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Figure 6. One of the two ground states of the competitive model, Eq. ([TJ, 
the Fr-AF region (Ji < 0, J2 < and \.h\ > ] J2I). 



Later on, it will be useful to consider the following decomposition of density free energy 
// of the non random model in the P phase (when necessary we will use the suffix / 
for stressing that the given quantity refers to the non random Ising model, as opposed 
to the random one) 



13 fi 



log2 + log(cosh(/3Ji)) + -log(cosh(/3J2)) +(^/(ti,t2). 



(12) 



The term ipj represents the non trivial part of the high temperature expansion of the 
density free energy, responsible for the critical behavior of the system. 

By using a recently proposed method [221 [301 [3T] , we will find the phase boundary 
of the random model by mapping it to the non random one. According to this mapping, 
to find the upper phase boundary of the random model we consider the transition 
equations ^ and ((HI) of the non random model and perform one of the two following 
substitutions 



ti = tanh(/3Ji) ^ t'-P = J dill (42) tanh Uji^}^ 
t2 = tanh(/3J2) ^ tf'> = J dii2 [J?)y) tanh h/^l 



(13) 



tanh(/3Ji 



(SG) 
1 



/ d^ll (42)tanh2(/3jii) 



t2 = tanh(/3 J2) ^ = / d^i2 ( J^l ) tanh^ ( /3 



(2)> 



(14) 



The substitution in Eqs. ([7^5]) provides the disordered ferro/antiferromagnetic 
transitions (P-F/AF), whereas the substitution (fH]) in Eqs. ([7]) provides the spin glass 
transition (P-SG). For a given couple of parameters ( Ji, J2), we find the corresponding 
critical inverse temperatures ff^^'^^^ or /3c^'^'', respectively. In the thermodynamic 



l{FIAF) 



and /3c 



(SG) 



limit, for given parameters (Ji, J2), between the two critical values /3c 
only one is stable and the system is critical at /3c given by 

/3c = min{/3f«),/3™}. (15) 

By performing the substitutions or (fH)) , the transition equations ^ and ^ 
take respectively the form 



4"^((4"^)' + 2C) + 2ir + l 



(16) 
(17) 
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where a = F for the substitution (jl3p and a ~ SG for the substitution (fill) . 

The multicritical point P-F-SG, if any, is solution of the system of equations 
obtained by taking the two versions of Eq. (|16|) with the two substitutions ([24]) and 
(|25p . Similarly for the multicritical point P-AF-SG from the two versions of Eq. (fT7|) . 

4-2. Critical behavior 

The transformations ([T^ and ([T^ can also be used to calculate the free energy and 
the correlation functions infinitely near the upper phase boundary above the P-F and 
the P-SG lines, respectively. To this aim, besides (y3/, we define also (p, that is, the 
non trivial part of the high temperature expansion of the density free energy / of the 
random model: 

^ /?/ = log2 + j dfii log(cosh(/3Ji)) + ljdfi2 log(cosh(/3J2)) + if- (18) 

According to [55], in the P phase, infinitely near the P-F and the P-SG lines, ip takes 
respectively the two following forms: 

^/(ii,t2)-¥'^^) =Pi{t^P,t^r), (19) 

^/(ti,t2)-^(^°)=i^.(tr\ir^). (20) 

Similar formulae hold for the correlation functions at zero external field. In other 
words, we get the exact critical behavior of the system in the upper critical line in the 
P region. On the other hand, as a general rule of the mapping, we immediately see that 
the critical exponents of the random system cannot be affected by the substitutions 
psp or (|14p . so that the critical behavior of the two systems in the P region is simply 
the same. Therefore, since the non random models defined over graphs of interest, 
like generalized Bethe lattices and Husimi tree (and many random graphs), have a 
mean field-like behavior with a second order phase transition, the same happens for 
the corresponding random versions. 

4-3. Condition for the mapping - the case of Bethe lattices and Husimi trees 

In |29| it has been shown that the upper critical boundary one obtains by using this 
mapping is exact whenever the underlying set of links is infinite dimensional in a broad 
sense as happens for instance on generalized tree-like structures having a finite number 
of loops per vertex. The condition for the mapping to be exact (infinite dimensional in 
the broad sense) is that ipi exists (see Sec. I) and that choosing randomly two arbitrary 
infinite long paths passing through a vertex, the probability p(/) that they overlap each 
other for I bonds goes to zero exponentially as / goes to infinity [H] . A path of length 
I is defined as a succession of I different bonds connecting vertices. In the framework 
of the high temperature expansion of the partition function, one sees that near the 
critical temperature, infinitely long closed paths correspond to the contributions of the 
free energy, while open paths and combinations of open and closed paths correspond 
to the contributions of the correlation functions. The mentioned probability for the 
overlapping of two paths is a feature that applies both to closed and open paths so 
that, for calculating p(Z) we are free to consider both open or closed paths. We can see 
easily that for our model p{l) — *■ exponentially as follows. Given I, let us consider 
two arbitrary paths of lengths li and I2 with li,l2 > Let us refer to these paths as 
first and second path, respectively. Since the paths are all statistically equivalent, we 
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Figure 7. Two paths starting from tlie vertex are drawn. Tiie bonds of the 
two paths are indicated as a\,a2, ■ ■ ■ and bi, i)2, • • •, respectively. The two paths 
overlap each other in correspondence of the bonds 02 = bg and 04 = 612. These 
bonds on the figure are slightly shifted for visual convenience. 

can choose for simplicity the first path to coincide with the straight fine indicated in 
Fig. El (focated in the most left part of the figure) . For the second path to overlap I 
times with I bonds of the first path, the second path must cross / vertices belonging to 
the first path. On the other hand, at any of such a crossing point, the second path has 
two over four possibilities (2/4) to make an overlap with a bond of the first path and, 
if the number of crossing points has to be infinite, as it must be for the requirement 
Z ^ 00, it is easy to see that only one of the above possibilities can be taken (choosing 
the bond going in the upward direction). Therefore, we see that p{l) < 1/2' and for I 
large we have p{l) ~ 1/4', i.e., an exponential decay. In general the same conclusion 
holds for any non trivial (that is with a vertex degree g > 2) generalized Bethc lattice 
and Husimi tree. Note that, in general, in such structures the number of closed paths 
passing through a vertex is infinite (see Sec. 1). 



5. Phase diagram of the random model 

In the two following subsections we apply the mapping to study the phase diagram of 
the disordered model introduced in Sec. 2.2. Note that, as manifested in Fig. [21 the 
P-F and the P-AF transitions are completely asymmetric one another and since they 
can occur only for positive and negative values of Ji , respectively, it will be sufficient 
to consider only the P-F and the P-SG transitions. Therefore, from now on we will 
study only the region with Ji > corresponding to Eqs. ^ and p6|) that here we 
rewrite 

4"'((4"^)'-2ti")) + l-24"' =0, a = F,SG. (21) 

We observe that Eq. ([2T|) (as well as Eq. (fT7|l ). with respect to arc quadratic 
and biquadratic under the transformations ([24|) and ([25]) , respectively, so that wc can 
solve this equation explicitly with respect to t^^\ This procedure turns out to be 
convenient for studying some asymptotes since we want to plot the critical lines on 
the plane (T/Ji, J2/J1). Let us put 



.(.) ^ l + ^-Vl + ^ + ^^.^ (22) 

X 

In terms of t^i \ the solution of Eq. (PT|) for 4'^'' can be expressed as 

ti^ = ^iA"^)- (23) 
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5.1. Triangular ± J model 



For the triangular ± J model the randomness is given through Eqs. ([5]) and 
plugged in the transformations and give respectively 



t\'> = {2p,-l)t, 

if' = (2p2 - l)t2 



and 



(SG) 
(SG) 



which, 



(24) 



(25) 



— 1-2- 

Note that in the last case there is no dependence on the values on pi and p2, so that 
in the triangular ± J model model the P-SG transition is unique. Note also that Eqs. 
(|24p are unchanged under the exchanges pi ^ \ — pi and Ji —Ji, and P2 ^ ^ — P2 
and J2 ^ ~^2, respectively. Therefore, it will be sufficient to restrict ourself to the 
region 

' Pi > 0.5, 
P2 > 0.5. 

From (|23p the solution for the P-SG line (in terms of t'"^'^^) is 



(26) 



n{SG) 



Jl 



tanh 



1/2 



while the solution for the P-F lines (in terms of t2 ) is 



rr^(J2 
Jl 



tanh 



2pi - 1 



(27) 



(28) 



As an argument based on the high temperature expansion of the free energy 
suggests, we see in Fig. |H]that the presence of loops, tuned by fixing Ji and increasing 
the amplitude of the parameter J2 , increases the spin glass critical temperature Tc 



(SG) 



In particular, the values of ri^'^-' for J; 



Tf°^(J2 



0) 



Jl 



tanh 



and J2 
V2 



00 arc given respectively by 



1.13459, 



(29) 



n(SG) 



(J2 



Jl 



tanh 



2-^/3 



= 1.74487. 



(30) 



As we shall sec below, for values of pi and p2 sufficiently high and far from 0.5 and 
for J2 > 0, a similar behavior holds for the disordered ferromagnetic transition as 
well. Note however that unlike Eq. fTT]) . for Eq. ((25)) exist regions without solution. 



fSG) (F) 

Finally, observe that, for Ji finite, the asymptotic values of Tc or Tc for J2 = 
and J2 = +00 correspond to the values of a pure Bethe lattice case and to a system with 
an infinitely strong loops effects, respectively. Note that in the latter case, if Ji 0, 
Tc and Tc also — ^ 0, the lattice being just an infinite set of non interacting dimers 
of spins: a gas of dimers. Wc point out that this last situation differs from the one in 
which J2 ^ ±00 and Ji remains finite: an arbitrary small but non zero coupling Ji 
make the structure with triangular loops non reducible to dimers. 
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5.1.1. The case pi = P2. Here we discuss an analysis carried out for the simpler case 
p = pi = P2. We have performed a study for the disordered P-F transition for several 
values of p. In Fig. [8] we report these lines. For p = 1 one recovers the non random 
ferromagnetic transition (Fig. [3]). In general for fixed parameters Ji and J2, as the 

(F) 

value of p is decreased from 1.00 toward 0.50, the critical temperature Tc decreases 
continuously. For a fixed value of p, as the value of J2 increases the P-F lines reach a 
horizontal asymptote whose value can be obtained explicitly from Eq. 



2 



tanh ^ 



A 



-1 



Ji 

where A = 2p — 1. When J2 decreases toward from Eq. (|28p we get 

tanh ' I — 
2A 



Ji 

Finally for negative values of J2 the P-F lines reach a final point given by Eq. 



(31) 



(32) 



Tf ^ -00) 
Ji 



tanh 



A 



(33) 



Note that, for fixed Ji, Eq. (|31|1 has no solution for values of p < 0.6909, whose 
asymptote corresponds to Tc /Ji = 0. Similarly, for fixed Ji and for J2 ^ —00 Eq. 
P3p has solution only for the pure ferromagnetic case p = 1. Note however that in 
both the cases the above temperatures Tc are lower than Tc so that, according 
to Eq. psp . the actual phase in these regions is SG and the above ri^'''s represent 
only unstable (not physically relevant) transitions. 

By using Eq. (fT5|) . in the plane (T/ Ji, J2/ Ji), as evident from Figs. [5] and 
[51 one finds the following scenario. Let us consider first the case J2 > 0. For 
p > 0.85 ± 0.01 for any value of the parameters Ji and J2 the transition can be 
only P-F. For p £ [0.81,0.85] ± 0.01 there is a range of values Ji and J2 where the 
stable transition is P-SG and a range where the stable transition is P-F. We note also 
that for p ^ 0.81 there can be more separated such ranges. Finally for p < 0.81 ± 0.01 
the transition can only be P-SG. Let us now consider the case J2 < 0. As is evident 
from Fig. whereas for J2/J1 > —1 for any value of p the system may have both 
a P-F or a P-SG transition, for J2/J1 < —1 the system has only a P-SG transition 
(similarly, for Ji < the P-AF transition may occur only if J2/J1 > 1)- Recall that, 
under the restriction of Eq. ([26|l . the averages of the couplings have the same sign 
of the parameters Ji and J2. For p > 0.5, the two completely different scenarios 
between the two cases J2 > and J2 < reflect the difference of the average of the 
frustration over the triangles, positive in the former and negative in the latter case. By 
looking back at the transformations ([M]) . we see that, concerning the disordered P-F 
transitions, the mapping consists in a simple renormalization of the couplings Ji and 
J2 so that the same argument used in Sec. 3 to describe the frustrated regions and, in 
particular, the dimers-triangles transition, applies again. On the other hand, unlike 
the non random case, for J2/J1 < —1 there exists always a stable and unique P-SG 
transition, signaling the onset of overlaps orderings. It is interesting to observe that a 
picture in terms of non interacting dimers can be found even in the glassy phase in the 

limit I J2I 3> I Ji| with Ji kept fixed. In such a limit in fact, the spins over a dimer must 

(2) (2) 

be exactly parallel or antiparallel, for Ji_^ > or Jx,y < 0, respectively, regardless 
of the values of the other couplings so that the average over the disorder of the 
local magnetization of any spin gives zero. 




Figure 8. Critical lines of the disordered transition P-F for several values of 
p = pi = P2 and the P-SG line. Note that, according to Eq. I|15|l . for any 
given p, there is a unique stable upper critical line obtained by taking the convex 
construction as shown in Fig. 1101 
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Figure 9. Enlargement of Fig. |8]around the SG line. 



In Fig. [TO]we report the detail of the phase diagram of the case p = 0.83 providing 
even the neighborhood of the F-SG crossover. 

Finally, we point out that, as it must be, for any choice of p, the multicritical 
point F-P-SG belongs to the Nishimori line. This fact can be immediately checked by 
substituting tanh(/3J) = 2p - 1 into Eqs. ([M]) and (^5]) . 

5.1.2. The general case. As shown in Fig. Illl for arbitrary values of the parameters 
pi and p2 the scenario of the possible transition lines becomes of course reacher. In 
the region J2 > we do not observe relevant qualitative differences with respect to 
the case pi = P2- We see that in general, for any fixed value of pi, by varying p2 one 
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T/ J: 




0.4 



,J2/Ji 



Figure 10. Phase diagram for the case Pi = P2 = 0.83 in the neighborhood of 
the multicritical point P-F-SG. The line drawn between the SG and the F phase 
is the Nishimori line which, as is known, passes inside the F phase. 



obtains a family of non intersecting critical lines margins at J2 = (as given by Eq. 
(|32p) and whose Tc in general becomes lower as p2 decreases. For fixed values of pi 
and p2, as the value of J2 increases, the P-F lines reach a horizontal asymptote whose 
value can be explicitly calculate from Eq. 



where A = 2pi ~ 1 and B ~ 2p2 — 1. From Fig. [TT]we observe also that two critical 
lines having two different values of pi, in general, may intersect each other. In the 
region J2 < we observe that if pi is sufficiently high, (starting from pi ^ 0.85) and 
P2 7^ Pi, as J2 decreases, the transition lines can have an inversion and reach another 
horizontal asymptote whose value is given by 



Note that the temperatures given by this equation are lower than the ones obtained in 
the opposite limit J2 +00. Note also that Eqs. and have no solution for 
values of pi and p2 such that A^B- 2AB - + 1 = and _ 2AB + 2A + l = 0, 
respectively. But analogously to the case previously seen (pi = P2), as evident from 
Fig. 1111 these regions of non existence occur only for values of pi and p2 quite lower 
than the threshold at which the only stable transition is P-SG. 

6. Conclusions 

The phase diagram of a competitive Ising model defined over the triangular Husimi 
tree has been analyzed, both in the non random and the random version (spin glass 
model). First, we have analyzed the non random case providing, in particular, an 
interesting physical characterization of the ground states in the frustrated regions 
in terms of a zero temperature phase transition between an ordered or antiordered 
state of connected triangles (F or AF) and a gas of non interacting two state dimers 
(also called a spin liquid in the literature). A natural question is whether or not this 
zero temperature phase transition represents the zero temperature boundary of a finite 
temperature phase transition. Second, we have introduced the disorder and studied the 
upper phase diagram and the corresponding critical behavior of the random case which, 




(34) 




(35) 
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Figure 11. Critical lines of the disordered transition P-F for several values of 
pi and p2 and the P-SG line. 

in particular, gives an afRrmativc answer to the above question. Our method relies 
on a recently proposed mapping valid for a wide family of graphs turning out to be 
infinite dimensional in a broad sense, as happens (but not only) in generalized tree-like 
graphs and in suitable structures where an infinite number of closed paths per vertex 
is also allowed. This second class includes, in particular, generalized Bethe lattices 
and Husimi trees, i.e., graphs over which non random models are exactly solvable. 
The mapping maps the spin glass model onto the corresponding non random one. For 
the two couplings J\ and Ji we have chosen the disorder giving rise to the triangular 
± J model. We can observe directly the role of the explicit frustration induced by the 
presence of the loops tuned by the amplitude and the sign of the parameter-coupling 
J2. As expected from an argument based on high temperature expansion and on 
frustration, the spin glass critical temperature is a nondecreasing function of IJ2I, 
whereas the disordered ferromagnetic critical temperatures is a nondecreasing function 
of J2. We observe the following scenario. In the non frustrated regions (that is the 
sectors I and III of the phase diagrams) , in the spin glass model we find an "ordinary" 
competition (like in the Sherrington Kirkpatrick model) between a P-SG and a P-F 
or P-AF transitions, in which the parameters J\ and J2 are simply renormalized by 
the mapping and the spin glass transition signals the onset of spin overlaps. In the 
frustrated regions instead (sectors II and IV of the phase diagrams), we find a new 
feature; in fact, even for | J2I > | Ji|, there exists always a phase transition which turns 
out to be a P-SG phase transition (the red line in the Figs. 8-11), whose glassy phase 
consists again of a gas of non interacting two state dimers. 

The analytical study carried out here concerns a non trivial model defined over 
a structure where loops are important. Other more complex and variegated cases, 
such as other Ising models with more than two couplings and Potts models, could be 
similarly analyzed. 

Even if, as stressed in Sec. I, in the Husimi tree, we do not have the existence of 
<p/, the density free energy in the thermodynamic limit, which is a necessary condition 
for the mapping, we think nevertheless that the application of the mapping also in 
this case provides important information of the random model. We stress in fact that 
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in the mapping no ad hoc approximation is introduced. 



Appendix A. 

A comment is in order concerning the model considered in the Ref. [13| (and partially 
in [m [H] ) , which is defined on the finite Husimi tree of Fig. O and the model given 
in Sec. 2.1, which is defined on the corresponding infinite Husimi tree (which can be 
seen as a Bethe lattice with loops). 

Whereas in the infinite Husimi tree one has to solve the model with the 
Hamiltonian ([T]) for vertices far from the boundary, and this is done by considering 
in all the calculation the lattice as infinite from the very beginning, in the finite 
Husimi tree one has to take into account also the boundary in the finite lattice and to 
perform the thermodynamic limit in the proper manner by calculating the correlation 
functions only from the density free energy after the thermodynamic limit has been 
taken. The relations between infinite and finite Husimi trees are on the same level of 
those between Bethe lattices and Cayley trees. It turns out that the physics, and in 
particular the free energy, described by these two versions of models is quite different 
[iQl mi [42] . Nevertheless, the equations defining the effective fields for the central 
spin used in both the approaches are the same implying that also the fixed points, 
and then the equations for the critical surfaces and the magnetization of the central 
spin, arc the same. Wc show this now by deriving the equations for the infinite Husimi 
tree following the so called Baxter approach [T], subsequently adapted to generalized 
Bethe lattices by many authors (see [10] and references therein) and to Husimi trees 

by [II]. 

Let us indicate with Z,j((To) the conditional partition function of the root spin 
(To (see Fig. [5]) for the n-th generation graph; that is the graph of Fig. [51 having n 
shells completed, having \Sn\ vertices, plus an external circular annulus C„, having 
the latter a number |C„| of vertices such that |S',i|/|C„| — *■ as n — *■ oo. We have 

Zr,=Y,Zn{<TQ), (A.l) 
Co 

where Z„ is the full partition function of the n-th generation graph. By looking at 
Fig. [21 it is easy to see that the conditional partition function Z„_|_i((7o) of the root 
spin at the n + 1-th generation graph is related to Z„((To) through 

Z„+i(ao) = e-''"°+'=^"°^"^+"^'+"^"^"^^n(<7i)Z„(a2), (A.2) 

where B = f3h, ki = /3Ji, k2 ~ f3J2, and cti and (72 arc the two Ising variables located 
at the vertices 1 and 2 respectively in Fig. [H It is convenient to introduce the effective 
field hn- 

Performing the sums over ui and U2 in Eq. (jA.2[) both for ctq = + and cto = — , and 
using the definition (|A.3[) we arrive at the following closed recursive equation for the 
effective field 

g2h„+i ^ 213 e +Ze +e 

g2fei+fc2 -f 2e~2fc2+2'ir. + e~2'=l + ''2+4'»n ■ V^'^l 

All the possible phase transitions (except the modulated ones) of the uniform model 
can be studied analyzing the fixed points h* of the above equation: 

„-2fci+fe2 I 9„-2fc2 + 2h* , „2fei+fe2+4h* 

2h* _ 2se j- i~ e I'A ^"l 

^ ~^ g2fci+fe2 2e~2'=2+2''* + e~2*i+'=2+4h* ■ l^-'^J 
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Of course one can use Eq. (|A.5P to write the recursive equation for the magnetization 
for the root spin (Tq- If one is interested in studying the individual magnetizations 
of further spins as for the spins cti and tT2, then one has to write the recursive 
equations for the conditional partition function Z„_|_i((To, cti, tT2) in terms of weighted 
sums over cri.i, . . . (72,2 of the conditional partition functions ^ri(f i, fi^i, (71^2) and 
■^n(c2, (72,1, 0-2,2) (see Fig. [5]) and defining the corresponding effective fields. In 
this way one can study possible modulated phases as done for example in [TUI for 
generalized Bethe lattices involving second and third-nearest neighbor interactions. 

In the model considered in this paper, the base spin erg is attached to one single 
triangle (see Fig. [5]). However one could be interested in studying cases in which the 
root spin has to be attached to more, say (7, triangles. It is immediate to recognize 
that in this case instead of Eq. (jA.Sp one has 

-2ki+k2 _|_ 2e-2fc2+2/i* _|_ g2fci+fc2+4/i* 



g2ki+k2 ^ 2e^2fe2+2/i* ^ g-2fei+fc2+4/i* 



The model considered in this paper corresponds to the case = 1 of Eq. (jA.sp which 
was also derived in [T3] for the finite Husimi tree. As already mentioned, despite 
the physics for the infinite and finite Husimi trees be different, the equations for the 
effective fields are the same, or in other words, the physics of the "central" spins, far 
from the boundary, are the same. Therefore we can, in particular, use the analysis of 
the Rcf. [13J for the phase transitions consisting in the study of the solutions of Eq. 
(jA.SP and their stability leading to the Eqs. (O and ©. 
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